The role of deformations in the bouncing droplet dynamics 
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Droplets bouncing on a vibrated liquid bath open ways to methods of manipulating droplets, 
creating double emulsion and performing pilot wave model experiments. Many periodic trajectories 
and resonant modes are observed when tuning the forcing parameters. This complex dynamics 
emphasizes the interplay between elastic energy storage and energy dissipation at each bounce. 

We propose to model droplets using a bouncing mass-spring-damper system that mimics a de- 
formable droplet bouncing on a non-deformable liquid bath. From the experimental measurements, 
we constructed bifurcation diagrams of the bouncing trajectories and challenged our bouncing spring 
model. The agreement between experiment and the spring model reveals that this model can be 
used to rationalize and predict a variety of bouncing droplets behaviors involving multi-periodicities. 

PACS numbers: 47.55.D-, 05 45.-a 



I. INTRODUCTION 

A droplet laid onto a vertically vibrated liquid bath 
can, under certain conditions, bounce periodically and in- 
definitely. Periodic trajectories are of crucial importance 
in a variety of bouncing droplets experiments. These 
droplets can present interactions at a distance [Iffi], 
deform following resonant modes [5], move using self- 
propulsion modes [SJ |S] or be used as experimental model 
for pilot wave systems [7H5]- These really interesting be- 
haviors rely on the periodicity of the bouncing droplet 
trajectory which can be multi-periodic or chaotic. The 
bouncing dynamics is strongly dependent on the relative 
bath/droplet deformation. However, there is a lack of un- 
derstanding on the droplet and the bath contributions to 
these bouncing systems. Here, we will answer the ques- 
tion of the droplet contribution by determining the role 
of the droplet deformation in the bouncing mechanisms. 

In this paper, we simplified the system by inhibiting 
the bath deformation and studying the role of the droplet 
deformation on the emergence of complex trajectories. 
We focus on small droplets bouncing on a highly vis- 
cous liquid bath that is vibrated vertically. The bounc- 
ing modes experimentally observed are summarized using 
bifurcation diagrams. We showed that a simple model 
like the bouncing ball is unable to explain the diversity 
of observed bouncing modes. We thus propose a model 
with the minimum physical ingredients required to cap- 
ture the main features observed experimentally. Indeed, 
we demonstrate that a bouncing droplet can be modeled 
as a bouncing Kelvin- Voigt material which is able to re- 
produce the trajectories of the droplet center of mass and 
their numerous periodicities when bouncing on a vibrat- 
ing bath. 

This simple model is able to rationalize and predict 
a variety of bouncing droplet trajectories on a non- 
deformable liquid bath that is oscillating vertically. This 
shall lead the path to future models that would incorpo- 
rate the bath deformation and be able to interpret spe- 
cific trajectories observed on a deformable bath such as 



the walking droplets [6 which are at the base of the pilot 
wave model experiments [THH] • 

We will start this paper by giving a global view of de- 
formations in the bouncing droplet problem. Then, we 
will describe the experimental setup and present an ex- 
periment of a droplet bouncing on a vibrating bath for 
a chosen frequency. The analysis of the the periodic- 
ity of the trajectories as a function of the control pa- 
rameter, namely the maximum acceleration of the plate, 
will be detailed. These results will be summarized using 
bifurcation diagrams. In the next section, we will pro- 
pose a minimal model based on a mass-spring-damper 
system to reproduce our experimental observations. Nu- 
merical simulations are performed to test the validity of 
the model. The numerical results are then presented in 
bifurcation diagrams which are compared to the experi- 
mental ones. Conclusions will close this paper. 



II. BACKGROUND 

In a series of recent works, bouncing modes have shown 
to exert particular features in bouncing droplet experi- 
ments. In both the pioneering work of Couder et al. |10j 
on the bouncing threshold and of Dorbolo et al. [5] on 
the various deformation modes, droplets are observed to 
bounce in phase with the vibrated bath. Couder et al. [6] 
discovered a self-propelled droplet, called 'walker', guided 
by the surface waves generated by the bounce on the sur- 
face of the bath. This peculiar behavior is made possi- 
ble when droplets double there bouncing period. Under 
these conditions, the droplets are able to trigger Fara- 
day surface waves and use them to propel themselves. 
By controlling the bouncing modes and consequently the 
surface waves produced by the bounce, Eddi et al. were 
capable of making different geometrical lattices of bounc- 
ing droplets [2] and making these bonded droplets move 
on the surface as ratchets [3]. 

Let us strive to have a global view of the problem. 
A liquid droplet of diameter D bounces on a liquid in- 
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FIG. 1: (Color online) (a) The asymptotic behavior of a bouncing droplet on a vibrating liquid bath is categorized and 
sketched as a function of the Ohnesorge parameters Ohd of the droplet and Ohb of the bath. These parameters represent the 
relative damping of the oscillations on the droplet and on the bath, (b) A selection of relevant works on bouncing droplets 
are plotted on the Ohnesorge diagram revealing the importance of the deformation of the droplet and/or of the bath in each 
experiment. 



terface oscillating sinusoidally at a frequency / and an 
amplitude A. The bath and the drop liquids are char- 
acterized by kinematic viscosities Vb and v&, densities pb 
and pa, and surface tensions Ofc and Od , respectively. In 
earlier works [11], we suggested that a relevant pa- 
rameter is the Ohnesorge number that compares viscous 
and capillary times. We then define Ohnesorge num- 
bers of the bath Ohb = Vb\Jp~bfob~D and of the droplet 
Ohd = Vd\/ pd/vd.D- These Ohnesorge numbers indicate 
how oscillations are damped by viscosity on the bath and 
on the droplet. Thus, these numbers give an indication 
on the relative importance of the bath and droplet defor- 
mations in experiment. On Fig. [T^t, we have represented 
a sketch illustrating the importance of these deforma- 
tions as a function of the Ohnesorge numbers Ohb and 
Ohd- Obviously, both parameters are not sufficient to 
describe univocally the behaviors of the droplets. This 
map on Fig. [lji can be seen as a projection of the muti- 
dimension phase diagram on a 2D map defined by Ohb 
and Ohd- 

On Fig. [JJj, we have situated a selection of bouncing 
droplet studies on the deformations map. The plotted 
areas were adjusted for a better visualization. The first 
model for the bouncing threshold has been developed by 
Couder et al. [JO] for highly viscous droplets on a highly 
viscous bath, both deformations of the bath and of the 
droplet are rapidly damped (top right corner of the map) . 
Using oil of viscosity 50 cSt, Protiere et al. investigated 
the interactions of droplets through the capillary waves 
they emitted and initiated fundamental studies on the 
walker's mechanism [H R)I [T2"rU4| . Eddi et al. performed 
pilot wave experiments using walkers made of 20 cSt sili- 



cone oil [H [15] . These series of experiments are based on 
small droplets (walkers) emitting waves on the liquid sur- 
face, these works are thus located on the left upper part 
of the map (cf. Fig. [IJd). Experiments on ratchet mo- 
tion, drifting rafts, Archimedean lattices and oscillations 
in crystal of droplets were performed with the same set 
of experimental conditions [IH3] • Fort et al. changed the 
viscosity of the droplet and the bath in order to change 
the length of the Faraday waves emitted by the walkers 
[5] , these experiments on the Landau levels are located in 
grey around the region described previously. During the 
preparation of this manuscript, we became aware that 
Bush et al. have submitted an analytical model which 
models the bouncing drop as a logarithmic spring on a 
low viscous bath [T6]. This work should be located near 
the walkers' region. Dorbolo et al. [5] focused more on 
the droplets deformation and used a highly viscous liquid 
bath to inhibit any deformation of the bath. We showed 
that a resonance phenomenon occurs between the bounc- 
ing droplets mode of deformation and the frequency of 
the vibrated bath, droplets deform as spherical harmon- 
ics which modes are selected by the forcing parameters. 
These deformation modes have a significant influence on 
the acceleration threshold for bouncing and this has been 
rationalized by an analytical model [11] . We also evi- 
denced that deformation of the droplet can be used to 
make droplets move horizontally on the surface of the 
bath [5], a non axisymmetric mode of deformation can 
be excited for specific forcing parameters which lead to 
self propulsion of the droplet on the bath. As well, these 
deformations can be used to create double emulsion in 
bouncing compound droplets [TT] ITS] . All these works 
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relying specifically on the droplet deformation are situ- 
ated in the lower right part of the map. 

In this work, we want to study the role of the droplet 
deformation in the bouncing droplet systems such as pi- 
lot wave experiment (walkers regions on Fig. [!}>). Thus, 
we decided to isolate the droplet deformation effect by 
positioning ourselves in the Ohnesorge diagram such as 
the droplet deformations are small and the deformation 
of the bath are negligible (highly viscous bath). This 
study should then be situated at the level of the Ohd of 
these walkers but at a relatively high Ohb (see blue disk 
on Fig. [1)3). That is where we will perform our study. 

III. EXPERIMENTAL SETUP 
A. Droplet generator 

We built a droplet dispenser to create submillimctric 
droplets of silicone oil (Dow Corning 200) of kinematic 
viscosity 20 cSt. The silicone oil surface tension is a = 
20.6 mN/m and density is p = 949kg/m 3 . On Fig. Ek, 
a small container, with a hole at the bottom, is closed 
by a piezoelectric chip at the top and is filled with the 
silicone oil. By injecting a short electric impulse (~ 5 ms) 
to the piezoelectric chip, a shock wave is released in the 
container and, through the hole, a drop is ejected which 
diameter is related to the size of the hole. In this case, 
the hole as a diameter of 600 /im and we are able to 
produce droplets which diameter is ranging from 650 /im 
to 750 /ira depending on the intensity of the impulse sent 
to the piezoelectric chip. For a impulse intensity value, 
the dispersion of the diameters of the droplets that are 
generated stays within 3%. With this technique, we are 
able to produce repeatedly droplets of silicone oil (and 
even water) of the same diameter in a range of 100 /im 
to a few millimeters by changing the hole diameter. For 
more information, see Ref 19J. 

In order to generate droplets of different diameters, 
one can also merge bouncing droplets on the bath. This 
can be made by pushing them against each other using 
the meniscus that is created by a stick dipped in the 
bath. For our study of the droplet trajectories, we chose 
to work with droplets of diameter 890 /im and made of 
20 cSt silicone oil which are typical droplets used for the 
20 cSt walking droplets |H] and pilot wave experiments 
[THS! [15]. Using the dispenser, we made two identical 
droplets of 700 /im on the vibrating bath. Then we force 
them to coalesce together and the resulting droplet has 
a diameter of D = 890 /mi. 

B. Vibrating bath 

The submillimctric droplets, of silicone oil 20 cSt, gen- 
erated by the droplet dispenser described above, are gen- 
tly laid on a highly viscous liquid bath made of 1000 
cSt silicone oil. This bath was vertically vibrated using 



an electromagnetic shaker GW-V55 following a sinusoidal 
motion of frequency / and amplitude A. We characterize 
the oscillation of the bath by the dimensionless forcing 
amplitude r = Air 2 Af 2 / g , where g is the acceleration of 
the gravity. 

An accelerometer (PCB-Piezotronics, 352C65), which 
delivers a tension proportional to the acceleration, was 
glued on the vibrating plate. With the latter linked to 
an oscilloscope, we measured the maximal dimensionless 
acceleration of the plate r. All the acceleration values in 
this manuscript are given with an error of 2 %. 

The frequency of the oscillating bath was arbitrarily 
chosen to be 50 Hz which is in the range of the forcing 
frequency commonly used in bouncing droplets experi- 
ments [J [21 H3 GO EJ E] and low enough to have a good 
time resolution while tracking with the high speed cam- 
era. 



C. Droplet tracking 

The bouncing droplets trajectories were tracked using 
the fast video camera N3 (1000 frames per second) from 
Integrated Device Technology, Inc. (IDT) (see Fig. [2)3). 
As fast video cameras demand a powerful lighting, we 
used several sets of 7LED-cluster from IDT. The posi- 
tion of the vibrating bath was tracked using a fixed wire 
glued to the plate which supported the bath. All the 
trackings, were made using the image analysis tracking 
tools provided by the software MotionProX from IDT. 
To have an estimation of the error on the measurements, 
the height of the bouncing droplet ranged from 0.5 to 3 
mm and the error on the tracking method was estimated 
to be about 0.05 mm. 



IV. EXPERIMENTAL RESULTS 

A. Bouncing modes 

For a fixed frequency of the oscillating bath and a forc- 
ing acceleration r higher then a threshold the droplet 
is observed to bounce periodically on the surface of the 
liquid bath as observed and described in previous publi- 
cations L3, 5, 10, 11 . At each bounce, the air film between 
the droplet and the surface of the bath is squeezed and 
can resist until the droplet bounces back into the air. 
This allows the air film to be renewed at each bounce. In 
some case, we can observe that droplets are permanently 
bouncing for acceleration lower than the gravity (I5 < 1). 
For our experimental set of parameters, droplets are sim- 
ply bouncing at the threshold, i.e. one bounce per oscil- 
lation of the bath, which deformation can be described 
as a spherical harmonic Y of degree 2 and mode (Y®) 
[11] . At each bounce the droplet loose energy in the 
droplet internal motion and in the lubrication film which 
are both eventually dissipated by viscosity. This loss at 
each bounce is compensated by the energy input coming 
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FIG. 2: (Color online) (a) Silicone oil droplets (20cSt) are 
made using a droplet dispenser which consists in a small con- 
tainer with a hole at the bottom and a piezoelectric chip at 
the top. A short electric impulse is injected to the piezoelec- 
tric chip which produces a shock wave in the container that 
ejects a droplet through the hole at the bottom, (b) These 
droplets were laid on the surface of a highly viscous silicone 
oil bath (1000 cSt) that is vibrated vertically using an elec- 
tromagnetic shaker. A high speed camera (1000 frames per 
second) recorded the motion and deformation of the drop from 
the side which were enhanced by a well positioned backlight. 



from the oscillation of the bath. This energy compensa- 
tion defines the threshold for bouncing. 

When r is further increased, above the acceleration 
threshold for bouncing Tj, the droplets experience com- 
plex trajectories that can be periodic over a few bath 
oscillations. Indeed, according to the oscillation phase, 
some extra energy can be used for a higher bounce and at 
the next bounce the energy given could be lower. How- 
ever, on average the energy loss at each bounce has to be 
balanced after a number p (an integer) of periods. The 
differences in bouncing heights are characterized by dif- 
ferent flight times between successive bounces which can 
be measured. We investigated the bouncing droplet tra- 
jectories as a function of the forcing parameters r of the 
oscillation of the bath, the forcing frequency / and the 
size of the drop D being kept constant. 

Using a drop of 890 fim on a bath oscillating verti- 
cally at / = 50 Hz, we recorded the trajectory from the 
side with a fast video camera at 1000 fps. Stable bounc- 
ing modes are sought by varying the initial conditions of 
the bounce for bath accelerations between = 0.9 and 



r = 5 by steps of approximately 0.1. In each movie, we 
recorded about 50 bath oscillations. The center of mass 
of the droplet and the bath positions were tracked over 
time. 

On Fig. |3j we present different characteristic periodical 
trajectories that corresponds to different forcing accelera- 
tions r. The position of the bath and the center of mass 
of the droplet are drawn in blue and red, respectively. 
The bath follows a sinusoidal motion. Trajectories of the 
center of mass of the droplet are vertically translated for 
better visualization. Note that the droplet is not bounc- 
ing instantaneously on the bath surface, it deforms itself 
when interacting with the bath that lasts about one third 
of a bath oscillation period. We define the droplet "con- 
tact" with the bath as the minimum in the trajectory of 
the center of mass of the droplet. 

Notable periodic bouncing trajectories can be observed 
on Fig. [3] with patterns that extend over a few bath oscil- 
lations and are periodically repeated. To characterize the 
different modes, we refer them by two parameters (p, q) 
according to the classical notification: the drop bounces 
q times while the bath oscillates p times. 

As we will see, the bifurcation diagram gives more in- 
formation on the bouncing modes than the (p,q) denom- 
ination. To construct the bifurcation diagram, we mea- 
sured two parameters characteristics of each bounce: the 
flight time Z\T m ; n and the height of the bounce h. We 
define AT m [ n as time intervals between two successive 
"contacts" of the droplet with the bath (that corresponds 
to two successive minima of the trajectory) and h as the 
height of successive bounces which is defined as the dif- 
ference of height in between a maximum and the previous 
minimum of the droplet trajectory. Both parameters are 
illustrated on Fig. [3^. The measurements of AT m \ n are 
normalized by the oscillation period of the bath T. We 
measured Z\T m ; n /T and ft as a function of the forcing _T 
on each of the experimental movies and reported these 
values as a function of r on Fig. [4^,-b. On these figures, 
red data points are related to the droplet trajectories 
observed as periodic while the grey data points are re- 
lated to the chaotic trajectories, no periodic motion was 
observed on these movies. The bouncing modes (p, q) 
are indicated above their corresponding red data points. 
These figures are called bifurcation diagrams as they de- 
pict the bouncing modes evolution as a function of the 
control parameter, F. 



B. Bifurcation diagrams 

The bifurcation diagrams on Fig. [4^,-b can be read as 
follow. A droplet permanently bounces on a bath oscil- 
lating at an acceleration r > Tf,. At the bouncing accel- 
eration threshold the droplet bounces simply on the 
bath following a mode (1,1). When r is increased, the 
bouncing droplet reaches a bifurcation when r = ^(2,2) j 
above which two successive bounces become uneven, cor- 
responding to the mode (2,2). After that, when r is 
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FIG. 3: (Color online) Experimental trajectories of a 20 cSt silicone oil droplet of diameter D = 890 /im bouncing on a 
bath oscillating at a frequency of 50 Hz for various accelerations _T. The bouncing mode (p, q) and the forcing acceleration r 
are indicated on each figure. The time interval ATmi n between two successive bounces and the bouncing heights h that are 
measured on each trajectory are illustrated on graph (e). 



increased further, a succession of bifurcations and modes 
(4,4), (3,2), (2,1) and (4,2) emerge and these modes can 
be interspersed with chaotic zones. 

To facilitate our understanding of the bouncing modes, 
we should analyze the bifurcation diagrams on Fig|4^i-b 
in parallel with the trajectories on Fig. [3j Indeed, these 
trajectories are represented on the diagram by one or 



several data points depending on the complexity of the 
trajectory. The first bouncing mode (1,1) is the mode 
which presents even bounces at each period of the oscil- 
lation of the bath. In consequence, it is characterized by 
the same time interval between droplet "contact" with 
the bath. However, two distinct heights can be observed 
which correspond to two different modes denoted by the 
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FIG. 4: (Color online) Trajectories characterization of the experiments and of the simulations (k = 0.072 N/m c = 
112 10" 6 kg/s) for a droplet of diameter 890 /im bouncing on a rigid liquid bath oscillating at 50 Hz for various forcing 
accelerations F. (a) Experimental measurements of the time intervals Z\T m i n , normalized by the oscillation period T of the 
bath, between two successive bounces; and (b) experimental measurements of the bouncing heights (cf. Fig. [3^). (c) Simula- 
tion measurements of the time intervals AT m i n /T; and (d) simulation measurements of the bouncing heights h. The different 
bouncing modes (p, q) are indicated on each diagram. 



same set of parameters (1,1) (see FigEja). The lower data 
points are characteristic of a mode (1,1) that we called 
(1,1)l, of the type presented on Fig. [3^, whereas the upper 
series of data points are of the type presented on Fig. [SJd 
called (1,1) h . The height of the (1,1) l type varies with 
r while the "contact" phase, observed on the trajectory 
plots, remains fixed. Conversely, the height of the (1,1) h 
type remains fixed while the "contact" phase adapts itself 
with increasing r. This latter case, is much more similar 



to a bouncing ball behavior [20 . As the acceleration in- 
creases, the droplet will continue to bounce in this (1,1) 
mode as long as the droplet can adjust the contact phase 
to compensate the loss of energy at each bounce. This 
locking of the phase is a well known feature in the bounc- 
ing ball phenomenon. However, the differentiation of the 
(1,1) mode into a high and a low mode is the consequence 
of the deformation of the droplet. This constitutes a main 
difference with the bouncing ball problem. 
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When r is increasing to r — 1.70 (see Fig. [3b), the 
droplet experiences a large bounce followed by a smaller 
one during two oscillations of the bath which mode is 
denoted by (2,2). When F is further increased, a wide 
variety of periodic modes can be observed. They are pre- 
sented on Fig. |3ji-h, the corresponding mode (p, q) and 
forcing acceleration r are written at the top of each fig- 
ure. These modes can be easily situated on the bifurca- 
tion diagrams on Fig|4^,-b. 

Another consequence of our deformable bouncing ob- 
ject is that the mode (2,1) is observed in two ranges of 
acceleration r » 3 and r « 4.3 (cf. Fig. |3k-b). Both 
modes described by the same parameters (2,1) are repre- 
sented on Fig. [3f and Fig. [3}i and occur for different r. 
They are characterized by a same time interval between 
the minima. However, the bouncing height of the mode 
occurring at higher _T increases with r as long as the 
contact phase remains constant. The other mode (2,1) 
is characterized by a constant height for an increasing 
r when the contact phase varies. This is qualitatively 
the same behaviors that is observed for the (1,1) h and 
(1,1)l modes. It is important to note that the mode 
(2,1) is known to be the characteristic mode of the walker 
droplets [?]. Even if walkers are observed on a low vis- 
cous (deformable) bath, the present work may suggest 
that several (2,1) modes may exist for the walker. 

From these observations, we can draw the following 
conclusions: 

• For a same set of forcing parameters, multiple tra- 
jectories can be stabilized. On Fig. [3^, and b, the 
mode (1,1)h and (1,1)l can be observed for the 
same parameter r, the selection of mode depend- 
ing on the initial conditions. The main differences 
between both modes reside in the bouncing heights 
and the contact phases which are different. 

• The classical notification (p, q) does not refer uni- 
vocally to one trajectory. Indeed, as for the mode 
(1,1), the mode (2,1) is observed for two different 
acceleration r = 2.70 (see Fig. M) and r = 4.35 
(see Fig. [3jr) and it is two different trajectories. 
This shows that the phase of the rebound is an ad- 
ditional relevant parameter. 

• According to the initial conditions, a trajectory at 
a fixed r can be either stable or chaotic. As two 
different stable modes, e.g. (1,1) h an d (1,1)l, can 
be selected depending on the initial condition, a 
chaotic or a stable mode can also be observed for 
a same set of parameters. Note that, under certain 
conditions, a droplet can alternatively pass from 
stable to chaotic and vice versa. 

We thus observed that the bouncing trajectories for a 
droplet on a vibrated bath is complex and modes are nu- 
merous. The bouncing ball model is too simple and is not 
able to reproduce the variety of trajectories that are ob- 
served [201 IS] ■ This is mainly due to the fact that in the 



bouncing droplet case, the droplet can deform. This ad- 
ditional degree of freedom must be taken into account as 
it influences the impact time which is not instantaneous. 
In order to take this effect into account, we propose to 
come back to a more basic system to determine the inter- 
relation between the impact and the deformation. For 
this purpose, we will start by a detailed observation of a 
droplet bouncing on a static bath. 



V. MODELLING 

A. Static liquid bath 

To gain insight on the bouncing mechanism of a 
droplet, we laid a droplet on the bath when it was at 
rest. We observed that the droplet bounces several times 
on the bath before coalescing with it. Using a high speed 
camera, we recorded an experiment, from the side, at 
2000 fps (frames per second). In this case a 740 /im 
droplet was considered. The reason of this particular size 
is that we study the behavior of a single droplet made by 
the dispenser. On Fig. [5] a succession of the droplet 
bounces is detailed. Based on the sequence of images, 
we constructed a spatio-temporal diagram according to 
the following method: the drop vertical centerlines of 
each image are juxtaposed. The time elapses from left to 
right. On the figure, one snapshot illustrates the droplet 
when it is squeezed at maximum during the impact. In 
this snapshot, the white vertical line indicates the pixel 
line represented on the spatio-temporal diagram. Be- 
fore impact, the drop is spherical and has a diameter 
D. Then, it deforms due to the impact. The maximal 
droplet horizontal extension is denoted by D + X. After 
that, the drop, which returns to a spherical shape, takes 
off, namely its center of mass goes up, without entering 
in contact with the underlying liquid. Residual oscilla- 
tions on the drop surface are then rapidly damped by 
viscosity. 

Physically the droplet has an initial kinetic energy 
that is converted in surface energy during the bounce. 
The surface energy is then restored and converted into 
kinetic energy of the center of mass of the droplet and 
into kinetic energy of the internal liquid motion. The 
latter is eventually dissipated by the viscosity while the 
kinetic energy of the center of mass is converted into 
gravitational potential energy which is maximal when the 
maximal height of rebound h is reached. 

Based on these observations, a model can be imagined 
to reproduce these bouncing modes. They are observed 
to be linked to the frequency and size of the droplets. A 
natural bouncing model would be the inelastic or semi- 
elastic bouncing ball which reproduces this kind of phase 
diagram. This model is characterized by an instanta- 
neous contact solely governed by the restitution coeffi- 
cient e of the ball. In our case, due to the finite contact 
time of the bouncing droplet, this model is too simple to 
capture the complex interaction in between the droplet 
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FIG. 5: (Color online) Spatio-temporal diagram of a 20 cSt 
silicone oil droplet of diameter D — 740 /im falling on a static 
highly viscous 1000 cSt silicone oil bath. Time elapses from 
left to right. The droplet experiences several bounces of 
heights hi, fi2 that are measured from the center of mass 
of the droplet when it is floating (white dotted line). The 
impact speeds of the first bounce is about 0.15 m/s. A snap- 
shot illustrate the shape of the droplet when its deformation 
D + X is maximal during the bounce. The center of mass is 
represented by green dots. The simulation of a falling mass- 
spring-damper system, from a height of hi — 1.1 mm, with 
parameters k = 0.072 N/m and c = 112 10 _6 kg/s, is rep- 
resented by blue and red lines. The central red line is the 
trajectory of the center of mass and both outer blue lines are 
the trajectories of the point masses mi and W2- 



and the oscillating bath. Let us remember that the term 
"contact" is improper because the droplet and the bath 
are never in contact. A droplet is considered to take off 
when the lubrication film is sufficiently thick for the force 
it exerts to be neglected. Gilet and Bush [531 [23] faced 
the same issue as they observed droplets bouncing on a 
vibrating soap film. They developed a theoretical model 
based on the second Newton law to predict trajectories. 
Depending on the shape of the soap film, they could de- 
duce the force acting on the droplet. Similarly, Eichwald 
et al. [24] investigated the phase diagrams of a solid 
bouncing ball on an elastic horizontal membrane whose 
frame oscillates. Note that in both works the bouncing 
motion is due to the deformation of the oscillating mem- 
brane or of the soap film and not from the deformation 
of the bouncing object. 

In a previous paper [5] , we demonstrated that a droplet 
behaves as a forced harmonic oscillator, the surface ten- 
sion providing the restoring force and viscosity the damp- 
ing. For Weber number at impact Wei < 1, the droplet 
stores energy in surface energy which scales as aX 2 [25] , 
where X is the length increment of the droplet at impact 
(cf. Fig.[5|. We can thus, naturally, propose a ID model 



for the droplet as a system of two masses linked by a 
spring which is characterized by a stiffness k in parallel to 
a dashpot with a damping coefficient c (cf . Fig. [6| . The 
behavior of the spring system is governed by the New- 
ton laws written for both masses. A one spring model 
without a dashpot has already been applied to infer a 
bouncing criteria and a restitution coefficient to inviscid 
droplets impacting an hydrophobic surface at high speed 
(corresponding to high deformation) [26 . On our side, 
droplets bounce on a vibrating bath with Weber number 
Wei < 1 or close to 1. Thus, we will study this sys- 
tem when the surface energy that is of the same order 
of magnitude as the droplet kinetic energy. In the next 
section, we will investigate the dynamics of the system 
bouncing on a static as well as on an oscillating plate. 
We started by numerically studying the behavior of a 
spring falling on a static plate. Then, we will compare 
the results to the experimental data (cf. Fig. [5]) in order 
to adjust the spring parameters k and c. Finally, we will 
study the bouncing modes of this mass-spring-damper 
system bouncing on an oscillating plate. We will then 
make a comparison with the experimental bifurcation di- 
agrams (cf. Fig. [4^,-b). We will show that, this model 
captures the main features that are experimentally ob- 
served which cannot be retrieved using the bouncing ball 
system governed by the single parameter, i.e. the coef- 
ficient of restitution. Indeed, in the case of deformable 
object, the coefficient of restitution is not constant and 
depends strongly on the dynamics of the bounce. 

B. Bouncing spring model 




FIG. 6: (Color online) The bouncing droplet is modeled by 
two masses mi,a linked by a spring of stiffness k in parallel to a 
dashpot with a damping coefficient c. The spring system is in 
contact with the plate oscillating sinuosidally at a frequency 
/ and amplitude A. As rri2 is in contact with the plate, a 
normal force N2 acts from the plate on ni2. In this case the 
spring is compressed (Ay < L), both masses feel outwards 
spring forces F s i and F S 2 in addition to the gravity forces F g i 
and F g 2- 

To reproduce the complex dynamics of the bouncing 
droplet, we model it by a Kelvin- Voigt material bouncing 
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on a oscillating rigid plate, as sketched on Fig. [6] The co- 
ordinate y is taken as the vertical one (positive upward) , 
Hi and y 2 are the coordinates of the upper mass mi and 
of the lower mass m 2 . The spring is characterized by a 
stiffness k and the dashpot by a damping coefficient c. 
At rest, the static length of the spring is L. 

On Fig. 6j we also represented the gravity forces F g i — 
—m\g anaF g2 — —m 2 g that act on both masses. The 
spring exerts also forces on each mass, they are F s \ = 
-k(Ay - L) and F s2 = k(Ay - L) with Ay = y x - y 2 . 
We add a viscous dissipation force which is proportional 
to the speed of elongation of the spring, with a viscous 
damping coefficient of proportionality c. When the lower 
mass m 2 is in contact with the plate, it experiences a 
normal force N 2 from the plate. 

The normal force N 2 depends in a complex manner on 
the elongation of the spring and on the bath accelera- 
tion which both vary during contact. Therefore we use 
molecular dynamics to compute this force N 2 when there 
is contact [37J|5S]. This normal force is null when there is 
no contact. Molecular dynamics simulation models nu- 
merically the contact between the plate and the mass 
m 2 as a damped spring interaction when it is in com- 
pression only. To model the fact that the plate is rigid, 
the mass of the plate m p is considered as a very large 
mass (m p 3> m 2 ) and the stiffness of the contact spring 
k c is chosen high enough to avoid any interaction with 
the spring that models the droplet (k c ^> k). Also, the 
damping coefficient modeling the contact is taken high 
enough to have an overdamped system modeling the con- 
tact. The value of these parameters, m p and k c , does not 
influence the results as long as they respect the above 
inequalities [29]. The contact dynamics is thus modeled 
as a very stiff over damped spring which gives a good 
approximation of the normal force exert by a rigid plate 
on a mass in contact. 

Considering the droplet model, Newton laws can be 
written for both masses mi and m 2 as 

m 1 y\ = -m x g - k{yi - y 2 - L) - c\yi - y 2 \, (1) 
m 2 y 2 = -m 2 g + k[y x - y 2 - L) + c\y x -y 2 \+ N 2 . (2) 

We choose to spread the mass of the droplet M on 
both masses, m x , 2 = M/2. The length of the spring L 
is set to 890 /im which is the typical size of our droplet. 
The spring constant k models the surface tension a 
0.02 N/m. This constant k determines the undamped 
frequency of the spring given by loq /2tt = \j2kjM. Us- 
ing the damping coefficient c, we can estimate the damp- 
ing ratio £ = — . c which compares the spring damping 
effect with the restoring effect. This damping ratio for 
the spring system is comparable to the Ohnesorge num- 
ber estimated for the droplet Ohd ~ 0.f5. The droplet 
spring system is thus underdamped and oscillates at a 
frequency u d = y/2k(l - £ 2 )/M. 

This model depends on two parameters only, k and c, 
which are the images of the parameters a d and re- 
spectively. We can fine tune them in order to reproduce 



the bouncing droplet trajectories as accurately as possi- 
ble. This is done by comparing the numerical simulation 
of the spring system laid on a static plate with the exper- 
imental results of a droplet laid on a static bath, which 
has been presented in Section fVA" 



C. Model parameters 

In the vibrating case, droplets bouncing periodically 
on a liquid bath oscillating at 50 Hz and F = 1 usually 
reach a maximum height of 0.5 mm which corresponds 
to an impact speed V* of about 0.10 m/s relatively to the 
system of reference of the laboratory. Moreover, they 
impact the moving bath whose speed is equal to Aui at 
maximum. As Aw is about 0.03 m/s for r = 1, the max- 
imal impact speed for those bouncing droplets is about 
0.13 m/s which correspond to a Weber number We w 0.6. 
This corresponds roughly to the 1st bounce shown on 
Fig. [5]for which Vi — 0.15 m/s. 

For velocities lower than \J a c ij pdF> ~ 0.17 m/s which 
correspond to We < 1, deformations are small during 
contact and the shape can be modeled by an oblate 
spheroid of major axis D + X [351 [3U]. When X < D, the 
increase in surface of the oblate spheroid is proportional 
to X 2 . We can write a simple energy balance between 
kinetic energy and surface energy, 



pD 6 V( ~ crA" 



(3) 



This latter equation justify the use of a spring to store 
energy during a bounce when We < 1. 

As the impact speeds of the bounces on the static 
bath (cf. Fig. [5]) are comparable to the relative impact 
speeds encountered for droplets bouncing on the oscillat- 
ing plate (cf. Fig. [3]), we adjusted the model parame- 
ters on the static case. On Fig. [5j the centre of mass of 
the droplet is measured and highlighted by green dots. 
Both spring-damper parameters k and c are adjusted in 
order to reproduce this succession of bounces as accu- 
rately as possible. In this perspective, we performed nu- 
merical simulations of the spring system laid from the 
height of T/i = hi + D/2 = 1.46 mm on a static plate 
(V, = y^ghl = 0.15 m/s and We, = 0.8). The pa- 
rameters values of the model are set in two steps, first 
the spring constant k is set using the deformation of the 
droplet and then the damping coeffiecient c is determined 
using the maximum height reached after impact. That is 
the spring constant k that mainly governs the maximal 
deformation of the droplet during the bounce. Firstly, we 
adjusted this parameter by setting c to and dropping 
the spring system on the static plate. By tuning the 
spring constant k to numerically reproduce the experi- 
mental deformations at impact, we set k to 0.072 N/m. 
Secondly, the damping coefficient c imposed mainly the 
height of the next rebound which height decreases when c 
increases. We set c to 112 10 _6 kg/s, which corresponds 
to £ = 0.33, in order to find h 2 the height of the next 
bounce (cf. Fig. [5]). We note that these two parameters 
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do not have independent effects. But this method was 
very satisfying to adjust the parameters in order to repro- 
duce the succession of bounces. We present the result of 
the numerical simulation on Fig. [5} The trajectory of the 
center of mass of the spring system is drawn in red while 
the trajectories of the masses mi and 777,2, assimilated to 
point masses, are drawn in blue. After comparing sim- 
ulations with experiments, we note that the trajectories 
of the center of mass are in good agreement. We observe 
that the floating drop (on the right hand side of Fig. [5]) 
has a lower center of mass than the prediction of the 
simulation. This could be due to the fact that the highly 
viscous bath is not completely rigid, so that the droplet 
create a small crater. Another difference is the additional 
oscillation of the spring system during the bounce which 
has, however, no influence on the next bounce. 



D. Numerical simulations 

Considering both calibrated parameters k and c for the 
model, we performed numerical simulations of the spring 
system bouncing on a 50 Hz oscillating plate. We have to 
note that the calibration was obtained for single droplets 
coming from the droplet generator while in our experi- 
ment, we tracked bouncing droplets that result from the 
coalescence of two bouncing droplets coming from the 
generator. Indeed the motivation is to study droplets 
which sizes are comparable to the Couder's walkers ob- 
served on a low viscous bath [31 01 [5] . Anyway, the values 
of c and k are not expected to vary with different sizes 
of the droplet. We adjusted the length of the spring L 
and the mass of the drop M according to the size of the 
droplet D = 890 /am that we used. This gives a damping 
ratio for the mass-spring-damper system £ = 0.25. 

We started by dropping the spring system, from a 
height of 1 mm, on a rigid plate oscillating at a frequency 
of 50 Hz and with an amplitude such that r = 1. The 
real time simulation was set to 11 seconds (w 550 oscilla- 
tions of the plate). We then increased r by steps of 0.1, 
without break between the simulations. In other words, 
the last state, the coordinates of the spring, its deforma- 
tion and the position of the plate, of one simulation is 
used as the initial conditions for the simulation with the 
next value of F. We repeated the process until reaching 
r = 5. After that, we decreased F by steps of 0.1 until 
reaching r — 0.5. We then processed the data of the 
center of mass of the spring system in the same way as 
we processed the data displayed selectively on Fig. [3] For 
each r value, we measured each bouncing height h and 
each time interval AT m - ln between minima of the center 
of mass trajectories. 

On Fig. |4j;, we reported the different intervals of time 
Z\T m j n between two successive bounces of the center of 
mass as a function of r. For each r value (increasing 
and decreasing phase), N measurements of Z\T m ; n were 
binned in classes of 0.1 ms of width (step of the simu- 
lation) and 0.05 of height (corresponding to the experi- 



mental error). We defined 77^ as the number of events in 
the bin i corresponding to time within the interval 0.1 i 
and 0.1 (i + 1) ms. On Fig. |4j:, the measurement recur- 
rences rii/N in each bin i are then represented in a grey 
level corresponding to (1 — rii/N) x 255 where a value of 
255 is white and a value of is black. A 100 % black bin 
means that all of the measurements corresponding to the 
acceleration r are located in the considered bin. Taking 
bins of 0.1 of width and 0.05 mm of height, we proceeded 
the same way with the measurements of h given by the 
numerical simulation. These are plotted on Fig. [4ji. We 
identified the different bouncing modes on both diagrams 
and indicated them on each figure. 

The agreement between the simulations and the ex- 
periments is remarkable. The first bifurcation from (1,1) 
to (2,2) is located at the same f(2,2) = 1-6 value. For 
F < F^2.2), the two different trajectories characterized 
by the same mode parameters (1,1), (1,1)h and (1,1) l, 
are observed in the simulations. The (1,1) degeneres- 
cence observed experimentally is thus reproduced by the 
model. The (3,2) mode is also observed and it is overly- 
ing the mode (2,1). In fact, the (3,2) mode was observed 
during the simulations when _T increased and the (2,1) 
mode was observed in the simulations during the decreas- 
ing phases. Thus, depending on the initial conditions one 
or the other mode can be observed. The bifurcation to 
the mode (4,2) is also observed for nearly the same accel- 
eration r in the experiment and simulation. The other 
mode (2,1) observed experimentally at T — 4.4 was de- 
tected when analyzing the trajectories of the simulations; 
however, it was not stable for more than several periods, 
which is why it does not appear clearly on the figure. 
The mass-spring-damper system model for the bouncing 
droplet bouncing is thus able to reproduce the key fea- 
tures observed in the experiments. These features would 
not have been evidenced by using a bouncing ball ap- 
proach. This demonstrates that the complexity and the 
variety of modes that are observed rely on the deforma- 
tion of the droplet and its finite contact time. 



VI. CONCLUSION AND PERSPECTIVES 

To have a global view of the various works on the 
bouncing droplets problem, we represented these on a 
2D map made with the Ohnesorge number of the bath 
Ohb and of the droplet Ohd- That map sorts the prob- 
lem as function of the relative importance of the defor- 
mation of the droplet and of the bath. In each case, 
the droplets trajectory can be complex depending on the 
forcing parameters and various bouncing modes can be 
observed. To highlight the influence of the droplet de- 
formation on the trajectory modes, we chose to focus 
on droplets, which deformations can be modeled as a 
quadratic spring, bouncing on a non-deformable liquid in- 
terface. This was made possible by considering droplet of 
diameter D — 890 /am, smaller than the capillary length 
A c ~ 1.5 mm, and when the bath liquid is highly viscous 
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(yi = 1000 cSt). We also chose to fix the frequency / to 
50 Hz which is in the range of frequency of recent works 
on the bouncing droplets [T| 121 IS] 151 HT1 IT?] . 

In these conditions, numerous bouncing modes with 
multi stable and chaotic regions can be observed when 
varying the acceleration r of the oscillating liquid bath. 
The variety of modes observed is richer than when a low 
viscous bath is used; in this later case the increase of r 
is limited by the appearance of the Faraday instability 
[2l 0]. We should note that the low viscous bath is a 
necessary condition to perform experiments with droplets 
interacting with the waves they emit when bouncing on 
the bath. To have further insight on the bouncing modes 
hierarchy, we made use of bifurcation diagrams in which 
we can have a global view of the modes succession, the 
multi stable and the chaotic regions. 

To reproduce the complexity and the variety of the 
bouncing modes, we proposed a model. An important 
feature of the bouncing droplets is that the droplet de- 
forms and stores energy in the surface energy before 
restoring it and taking off, this leads to a finite contact 
time. That is why, the bouncing ball model cannot cap- 
ture the behaviors that we observe experimentally. We 
thus proposed a model that can capture the deformation 
of the droplet and consequently the finite contact time 
with bath. It consists in two masses linked by a spring 
of stiffness k in parallel with a dashpot of damping coef- 
ficient c. These two parameters could be adjusted dur- 
ing experiments which were performed on a static bath 
before extending the comparison to an oscillating bath. 
We thus simulated the spring system on an oscillating 
rigid plate and found the bifurcation diagrams that we 
observed experimentally. The agreement between the ex- 
periments and the simulations is good. This model can 
even retrieve the multi stable state like the (1,1) h mode 



and the (1,1) l mode degenerescence. It can also cap- 
ture the two different modes characterized by the same 
parameter (2,1) which are observed for different accelera- 
tion r. This shows that the mass-spring-damper system 
(Kelvin- Voigt material) is an excellent model for bounc- 
ing droplets. 

Regarding the mode (2,1), we can imagine that dif- 
ferent (2,1) modes may also occur in the case of the 
'walkers' discovered by Couder et al. [BJ. Enlightened by 
this work, we should take caution and investigate further 
which of these 2 modes (2,1) is able to trigger the Faraday 
waves and make the droplet walk on the liquid surface. 
In the future, it would be very interesting to test the de- 
pendence of this model with the frequency. Indeed the 
deformation of the droplet and the mass-sprin g-damper 
system occurs in a characteristic time t c ^ ~ yl pdD 3 /a b 
and t c , s ~ y/M/k, respectively. We should thus observe 
and show a dependence of the bouncing modes on the 
frequency. Later, the model shall be extended to low vis- 
cous bath. The reason is that the coupling between the 
bath and the droplet mode are the origin of fascinating 
behaviors of the bouncing droplet like the walkers [6J , the 
ratchets [3] and the droplet crystals [2J. 
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